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ABSTRACT. This paper examines altitudinal changes in butterfly communities in 
the Sierra de Javalambre of central Spain. From May to September 1991 we sampled the 
butterfly fauna every 15 days at 10 stations located between 1100 and 2000 m. We 
obtained a total of 2,123 individuals of 101 different butterfly species. Abundance and 
species richness were highest at low elevations and declined with increasing altitude, 
while the converse held for habitat breadth. The chorological index was highest at low 
elevations. Changes in the butterfly communities along the altitudinal gradient of the 
Sierra de Javalambre are caused by harshness of envirommental conditions, changes in 
the vegetation, and presumably, resource impoverishement. High elevations do not seem 
to ’select’ for a endemic fauna of butterflies. The communities in the lowest places are 
composed of rare and localized species, while high elevation communities are less original 
in faunistic composition since they are composed of euryecious and widespread species 
in this area. 

Additional key words: altitudinal changes, habitat breadth, resource poverty, species 

richness. 



Analyses of altitudinal changes in diversity, abundance, and species 
composition of biotas can provide important information on such phe- 
nomena as those aspects of the enviromment limiting the distribution 
of organisms, factors influencing the structure of communities, and 
biogeographical patterns. These problems have been investigated by 
ecologists for the past 40 years. Several studies (Hagvar 1976, Claridge 
& Singhrao 1978, Hebert 1980, Ichijo et al. 1982) have concluded that 
a decrease in species richness with elevation is typical of many groups 
of animals, including insects, with the exception of bees (Gauld 1987) 
and tropical psocids (Turner & Broadhead 1974). At least four causes 
have been suggested for this decrease: reduced habitat area at high 
elevations, reduced resource diversity, increasingly unfavorable envi- 
ronments, and/or reduced primary productivity at higher altitudes 
(Lawton et al. 1987). In addition to these, other processes (competition, 
predation, evolutionary time, etc.) also may influence species richness 
(Lawton et al. 1987). 

Other studies (Janzen 1973, Janzen et al. 1976), however, have con- 
cluded that species richness peaks at middle elevations, rather than at 
low ones. Many distinct processes have been proposed to explain mid- 
elevational peaks: disturbances caused by human activities in low ele- 
vations (Wolda 1987), increasingly unfavorable environments at both 
high and low altitudes (Gagne 1979, Randall 1982), and even the sam- 
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pling regimes (short-term vs. long term sampling regimes) (Wolda 1987, 
McCoy 1990). 

The purpose of this paper is to present the results of a study of 
changes in the structure of butterfly communities along an altitudinal 
gradient in a Mediterranean mountain. There is a dearth of such studies 
in these latitudes (see Claridge & Singhrao 1978). 

Materials and Methods 

The study was conducted in the Sierra de Javalambre, Sistema Iberico, 
Teruel province, central-eastern Spain (Fig. 1). The Sierra consists large- 
ly of carbonated Jurassic rocks. The climate is typical of mediterranean 
mountains, with a wide daily and seasonal variation in temperature, 
with hot summers and cold winters. Rainfall is low and seasonal, varying 
anually from 400 to over 500 mm. These mountains thus provide a 
remarkable opportunity for the study of altitudinal distributions and 
zonation of organisms. 

The lower slopes are heavily cultivated with various fruits and veg- 
etables. Elsewhere, forests dominated by holm oak {Quercus rotundi- 
folia Lam.) (Fagaceae) and Spanish juniper {Juniperus thurifera Lin- 
naeus) (Cupressaceae), and their successional stages dominate the land- 
scape. At about 1300 m the forest is dominated by lusitanian oak {Quer- 
cus faginea Lam.) (Fagaceae), extremely perturbed and replaced in 
many cases by bushes of Ligustrum vulgar e Linnaeus (Oleaceae), Pru- 
nus spinosa Linnaeus (Rosaceae), Berberis vulgaris Linnaeus (Berber- 
idaceae) and Rosa canina Linnaeus (Rosaceae). 

Spanish juniper and black pine (Pinus nigra Arnold) (Pinaceae) form 
mixed forests at 1400-1600 m, dominated by seminatural black pines. 
This forest is replaced above 1600 m by Scot pine forests (Pinus syl- 
vestris Linnaeus) (Pinaceae). Just below the summit, at about 2000 m, 
there is a meadow with Erodium celtibericum Pau (Geraniaceae). A 
general account of the natural vegetation of the area is given by Peinado 
and Martinez-Parras (1985, 1987) (also see Fig. 2). 

Methods 

The study was conducted along an altitudinal transect between the 
village of Camarena de la Sierra and the peak of the highest mountain 
(Javalambre, 2020 m). Ten stations separated by intervals of approxi- 
mately 100 m of elevation (from 1100 to 2000 m) were established in 
the above described habitats. 

Every two weeks from May to September 1991 we sampled the 
butterfly fauna at these 10 stations. Samples were taken on sunny days 
between 1000 and 1600 h. The sampling scheme was based in sampling 
subunits (40 min) of collecting effort per site per sample. We obtained 
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Fig. 1. Map showing the study area. Transect with the numbers of sampling sites are 
also represented. 



a total of 2,123 individuals of 101 different species (Appendix A). From 
these data, the following variables were calculated: 

— Number of species (S) equals the total number of species recorded on each site. 

— Abundance (Ab) equals the total number of individuals. 



Volume 49, Number 3 



195 





Fig. 2. Schematic representation of both potential (A) and present (B) vegetation 
along the elevational transect. Numbers 1-10 refer to sampling sites, number 1 being at 
the highest elevation. 
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Table 1. Values of the 7 variables used in the characterization of butterfly com- 
munities on each sampling site, number 1 being at the highest elevation. 
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— Dominance (D) was obtained from the McNaughton & Wolff’s index (1970), expresed 
as D=Y,,/Y, where 2 is the sum of individuals of the two most abundant species, 
while Y is the total number of individuals. 

— Diversity (H') was obtained from the Shannon & Weaver’s index (1963), expresed as 
H'= -S Pi log2 Pi where Pi is the proportion of the ith species in the total sample 

— Equitability (J') expresed as: J'=HVH'max. 

— Mean chorological index (Cl). Calculated for each site by means of the chorological 
index of each species (data from Kudrna 1986). As Kudrna pointed out, the chorological 
index allows an evaluation of the biogeographic arrangement of all European butterfly 
species from a conservation point of view, as well as an evaluation and comparison of 
habitats (localities) based solely upon the composition of their butterfly fauna. The 
values of Cl result from the sum of the numerical values of range size, range composition 
and range affinity (“sensu” Kurdna 1986). Values for range size vary from 1 (species 
widespread across Europe) to 5 (species confined to small areas, such as islands, mountain 
ranges, or single sites in Europe). Range composition evaluates the continuity of the 
distribution; in other words, the ability of individuals of one population to reach other 
populations. It ranges from 1 (continuous, or nearly so, distribution over most of the 
European range of the species) to 5 (widely isolated single populations, small groups 
of populations, and small stocks of very restricted range). Lastly, range affinity syn- 
thesizes the relationship between the species’ European distribution and its world range 
as a supplementary indicator of the relative importance of these populations for the 
overall survival of the species. Its value ranges from 1 (extra-European species as defined 
by Kurdna 1986) to 4 (species endemic to Europe). Thus, the chorological index varies 
from 4 (most widespread species) to 14 (endemic European species restricted to very 
small territories). 

— Mean habitat breadth (HB). Obtained from the habitat breadth of each species after 
Simpson’s (1949) formula, expressed as HB= l/2p^, were Pj is the proportion of in- 
dividuals on each site. 

The values of all these variables are summarized in Table 1. 



Results 

Table 2 shows the results of the regression analysis between altitude 
and variables of community structure. Significant negative relationships 
were found between altitude and both species number and abundance 
(i.e., species number and abundance decrease with altitude), whereas 
mean habitat breadth increased with elevation (Fig. 3). However, no 
correlations between altitude and dominance, species diversity or equit- 
ability were found. This is due to the fact that number of species (species 
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Table 2. Results of the regression analysis between altitude and the 7 variables 
considered. Cl* represent results of the regression after the localities 8 and 10 were 
removed (see text for explanations). 
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* P < 0.05. = Coefficient of determination (in %). 



richness) and equitability (evennes in abundance) are the two com- 
ponents defining H'. Diversity increases as species are added, as well 
as when the species abundances are evenly distributed. In diverse sit- 
uations, single species do not dominate; in contrast, where one or two 
species are much more abundant than the rest, there is low diversity 
(see Price 1984). Thus, the absence of a relationship between elevation 
and equitability is the result of the absence of a relationship between 
elevation and butterfly species diversity as expressed by H'. 

As regards the mean chorological index, no significant relationships 
were found when all 10 sampling sites were included in the analysis. 
However, when localities 8 and 10 (both located close to human-per- 
turbed habitats) were removed from the analysis, a highly significant 
negative relationship between elevation and mean chorological index 
was found (Fig. 4). 

To evaluate the effect of the resource impoverishment (i.e., decrease 
in the abundance of suitable host plants) on the variation in species 
numbers along the gradient, we proceeded as follows. We assigned each 
butterfly species to one of four groups based on their larval host plants: 
1) plants of the class Dillenidae (families Brassicaceae, Resedaceae, 
Primulaceae, Malvaceae); 2) plants of the class Rosidae (families Ro- 
saceae, Fabaceae, Rutaceae, Umbelliferae); 3) plants of the class As- 
teridae (families Caprifoliaceae, Valerianaceae, Dipsacaceae, Boragi- 
naceae, Lamiaceae, Escrofulariaceae); and 4) plants of the class Lillidae 
(family Poaceae) (Appendix B). Figure 5 shows that in the lowest el- 
evations the percentages of the four different butterfly groups are close 
to those expected (high values of equitability), while at the highest 
altitudes some groups are proportionally better represented than ex- 
pected (low values of equitability). 

Discussion 

Studies on the distribution of insects along elevational gradients have 
yielded differing results (see McCoy 1990 for a review). Recent long- 
term sampling studies have concluded that previously identified mid- 
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Fig. 3. Relationships between altitude and species number (A), abundance (B), and 
mean habitat breadth (C). 
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Fig. 4. Relationships between altitude and the mean chorological index. Dashed line 
obtained from all 10 sampling sites. Bold line obtained after sampling sites 8 and 10 were 
removed (see text for explanations). 



elevational peaks in species richness may have resulted from the short- 
term sampling regimes that were employed (Wolda 1987). However, 
mid-elevational peaks have been documented for organisms that should 
be less affected by factors such as sampling regime and disturbance, 
for example, birds of the Paramos (Vuilleumier & Ewert 1978, Terborgh 
1977). 

Our data indicate that both richness and abundance of butterflies 
decrease with altitude, although the greatest number of species was 
found at 1500 m. This may be interpreted as a result of the interaction 
of two different variables-habitat and disturbance. The lowest places 
of Javalambre have been human-managed historically for agricultural 
uses, thus the climax vegetation has been replaced by its serai stages, 
which support lower butterfly species richness (Baz 1986). On the other 
hand, habitats at mid-elevations are deciduous woodlands {Quercus 
faginea) that contain the richest butterfly faunas, at least in the Iberian 
peninsula (Baz 1987, Viejo et al. 1989). If the habitat at lower elevations 
were in a natural stage, one would expect to found more butterfly 
species there. The continuous decrease in the number of species with 
increasing elevation may be caused by the harshness of envirommental 
conditions and area reduction, but also as a consequence of a reduction 
in resource diversity (Lawton et al. 1987). In the absence of more 
detailed inventories of plant species. Fig. 5 may serve as an example 
to illustrate this point, showing that at the highest elevations a trophic 
group of species dominates community composition (ca. 55% of all 
species). 

An interesting pattern found in this work is an extension of Rapoport’s 
latitudinal rule to altitude (i.e., the tendency for latitudinal ranges to 
become smaller with decreasing latitude [Rapoport 1982, Stevens 1989]) 
which has been reviewed recently by Stevens (1992). Mean habitat 
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Fig. 5. Altitudinal changes in the relative importance of the 4 butterfly groups in 
relation to their caterpillar host plants. The altitudinal changes of equitability of the 
butterfly assemblage are also represented. Sampling site number 1 is at the highest 
elevation. 



breadth increases with altitude (Fig. 3), hence the butterfly communities 
are composed primarily of common and euryecious species in the high- 
est elevations (see Baz 1987 and Thomas & Mallorie 1985a who found 
the same pattern in other Mediterranean mountains). Stevens (1989) 
indicates that the latitudinal RapoporFs rule arises as a result of the 
narrowing range of climatic conditions the individuals experience with 
decreasing latitude. Since non-migratory individuals must be physio- 
logically or behaviorally capable of tolerating the full range of condi- 
tions the seasonal changes impose on them, natural selection has favored 
broad tolerances at high latitudes. The consequence is that species from 
high latitudes have large latitudinal ranges because each individual of 
the species must have broad climatic tolerances just to survive at any 
latitude location. 

At the other extreme, the breadth of climatic conditions each indi- 
vidual of a tropical species experiences is so narrow that there is no 
penalty for possessing narrow climatic tolerances. As Stevens (1992) 
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pointed out, if this explanation for Rapoport’s latitudinal rule is correct, 
a pattern conforming to the altitudinal Rapoport’s rule should be found, 
since the breadth of climatic conditions experienced by mountain res- 
idents does increase with increasing elevation. Certainly, our data sup- 
port this idea. 

Lastly, an interesting pattern in relation to the value of mountains 
for conservation purposes has been observed. Figure 5 shows that after 
man-managed localities (i.e., 8 and 10) are eliminated from the analyses, 
an inverse correlation between the mean chorological index and altitude 
exists, so that communities in the piedmont are composed of more rare 
or localized species (in a European context) than those at the highest 
elevations which probably is due to the peculiarities of basal habitats 
(Baz 1987, Thomas & Mallorie 1985b). In Iberia, mountains chains have 
acted as corridors for butterfly species of predominantly European 
origin (Martin & Gurrea 1990), and have facilitated the expansion in 
range of many European species. As a consequence, low altitudes seem 
to have inhibited dispersal of banal European butterflies (Martin & 
Gurrea 1990) and nowadays harbor a more original Mediterranean 
butterfly fauna. However, more detailed studies in undisturbed habitats 
within altitudinal ranges 8 and 10 of our study are neccessary to prove 
this hypothesis. 
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123456789 10 Total Cl HB 

Papilionidae 

Iphiclides feisthamelii (Duponchel, 1832) 1 21421 118 4.48 

Papilio machaon Linnaeus, 1758 2 1 1 4 5 2.66 
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Hyponephele lupina (Costa, 1836) 1 1 1222 99 5. 4 
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Kanetisa circe (Fabricius, 1775) 11 11237 16 9 3.87 
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Appendix B. Butterfly species included in each group of caterpillar host plant. Sources 
from Takhtajan (1980), Gomez-Bustillo and Fernandez- Rubio (1974) and Higgins and 
Riley (1980). 



Dillenidae 


Rosidae 


Asteridae 


Lillidae 


A. cardamines 


1. feisthamelii 


P. panoptes 


P. hispana 


A. euphenoides 


P. machaon 


A. urticae 


A. arethusa 


E. ausonia 


P. apollo 


E. aurinia 


C. briseis 


P. brassicae 


A. Gratae gi 


E. desfontainii 


C. arcania 


P. rapae 


C. alfacariensis 


L. reduct a 


C. dorus 


P. daplidice 


C. crocea 


M. cinxia 


E. epistygne 


L. alciphon 


G. cleopatra 


M. didyma 


E. triaria 


L. phaleas 


G. rhamni 


M. parthenoides 


E. zapateri 


A. adippe 


L. sinapis 


M. deione 


H. fidia 


A. aglaja 


A. allous 


M. athalia 


H. hermione 


A. niobe 


A. cramera 


M. phoebe 


H. semele 


A. pandora 


C. rubi 


V. cardui 


H. statilinus 


A. paphia 


C. argiolus 


C. boeticus 


H. lupina 


I. lathonia 


C. osiris 


S. proto 


H. lycaon 


C. alceae 


C. semiargus 




K. circe 


P. alveus 


G. alexis 




L. megera 


P. carthami 


G. melanops 




M. jurtina 




L. boeticus 




M. ines 




L. idas 




M. lachesis 




N. acaciae 




M. russiae 




N. spini 




P. aegeria 




P. argus 




P. bathseba 




P. pylaon 




P. cecilia 




P. amandus 




P. tithonus 




P. bellargus 




S. acteae 




P. damon 




P. onopordi 




P. dorylas 




T. acteon 




P. fabressei 
P. icarus 
P. nivescens 




T. flavus 




P. ripartii 
P. thersites 








S. orion 
S. pirithous 
L. celtis 
B. hecate 
E. tages 
H. comma 
P. armoricanus 
P. cirsii 
P. serratulae 
S. sertorius 







